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Abstract. By resorting to the thick-chain model we discuss how the stretching response 
of a polymer is influenced by the self-avoidance entailed by its finite thickness. The 
characterization of the force versus extension curve for a thick chain is carried out through 
extensive stochastic simulations. The computational results are captured by an analytic 
expression that is used to fit experimental stretching measurements carried out on DNA 
and single-stranded RNA (poly-U) in various solutions. This strategy allows us to infer 
the apparent diameter of two biologically-relevant polyelectrolytes, namely DNA and poly- 
U, for different ionic strengths. Due to the very different degree of flexibility of the two 
molecules, the results provide insight into how the apparent diameter is influenced by the 
interplay between the (solution-dependent) Debye screening length and the polymers' "bare" 
thickness. For DNA, the electrostatic contribution to the effective radius, A, is found to be 
about 5 times larger than the Debye screening length, consistently with previous theoretical 
predictions for highly-charged stiff rods. For the more flexible poly-U chains the electrostatic 
contribution to A is found to be significantly smaller than the Debye screening length. 
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Introduction 

In recent years, the remarkable advancement of single-molecule manipulation techniques 
has made possible to characterise with great accuracy how various biopolymers respond 
to mechanical stretching QI2I3I1I21IID0I3I3GI3III1II2II3I14). The wealth of 
collected experimental data have constituted and still represent an invaluable and challenging 
benchmark for models of polymers' elasticity ||T5][Tl|[T71[IlliniliniEIlE2llIlE31- The 
interpretation of single-molecule stretching experiments often relies on one-dimensional non- 
self-avoiding models of polymers. It is physically appealing that the schematic nature of 
such descriptions often conjugates with the capability of reproducing well experimental 
measurements. Two notable instances are represented by the freely-jointed chain 1251 1151 
and the worm-like chain models which, in their original or extended forms, constitute the 
most commonly-used theoretical frameworks for biopolymers' stretching 1 26l 1171 l27l . Both 
models are endowed with a parameter, the Kuhn length or the persistence length, that provides 
a phenomenological measure of the polymer stiffness and that is obtained by fitting the 
experimental data. It is important to notice, however, that it is possible to go beyond this 
phenomenological approach and connect the persistence length to the fundamental structural 
properties of a polymer. A strong indication of the feasibility of this scheme is provided 
by the fact that, for a large number of biopolymers, the observed persistence length shows 
an approximate quartic dependence on the polymer diameter, as predicted for "ideal" elastic 
rods 1 131 . 

From this perspective it appears natural to investigate in detail the connection between 
structural properties and stretching response of biopolymers. We have recently pursued this 
objective by modelling explicitly the intrinsic thickness of a homo-polymer (treated as a tube 
with uniform cross-section) and characterising the stretching response 1211 . The theoretical 
and numerical results were employed in an appealing chemico-physical framework where the 
diameter of a biopolymer was not probed directly but inferred through the mere knowledge 
of the stretching response. The adopted thick-chain model |28 29 30 1, briefly outlined in the 
next section, was used to fit stretching measurements obtained for a variety of biopolymers: 
DNA 1 5 1, the PEVK-domain of the titin protein |7 8 1 and cellulose 13 1111 01 - For uncharged 
polymers, such as titin and cellulose, the effective diameter recovered from fitting the force- 
extension curves were very consistent with the stereochemical ones, thereby validating the 
thick-chain model approach |21 1. Even more interesting is the case of polymers possessing a 
substantial linear charge density, such as DNA and RNA which will be the focus of the present 
study. The properties of polyelectrolytes, in fact, depend very strongly on the electrostatic 
screening provided by the ions present in solution. The influence of the electrostatic screening 
on the behaviour of polyelectrolytes has been extensively investigated both experimentally 
and theoretically 11321 1331 l34l 1351 1361 [37 38 1. From the latter perspective, it is customary to 
introduce apparent (or effective) physico-chemical parameters to describe the properties of a 
polyelectrolyte in a given solution with reference to the uncharged polymer case. For example, 
in the context of elasticity, one introduces an effective (solution-dependent) bending rigidity 
to account for the additional electrostatic contribution to the "bare" persistence length of the 
hypothetically-neutral polymer |37 38 1. Also, in the context of colloidal dispersions of stiff 
polyelectrolytes, one can describe the polymer as uncharged cylinders and resort to the theory 
for second virial coefficient to derive its solution-dependent effective diameterl 34l 1391 l40l . 
For DNA in solutions of low ionic strength, both the effective persistence length and effective 
diameter can exceed by several factors the bare ones. So far, these effective DNA properties 
have been typically probed by distinct methodologies. For example stretching experiments 
were employed to establish the dependence of the the persistence length on ionic strength 
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1411 [&j while measurements of second virial coefficients, knotting probabilities or braiding 
properties were used for the effective diameter| 42 40 43 44 45)). 

The thick-chain framework is used to obtain, starting from stretching measurements 
data, the effective diameter of a polyelectrolyte and to further relate it to its the effective 
persistence length. Besides the implications in the general context of polyelectrolytes the 
proposed method can be used to establish the effective structural parameters to be used in 
coarse-grained studies of looping, knotting and packaging of biopolymers I46ll47l. 



The thick chain model 



To characterise the stretching response of a polymer with finite thickness we shall view the 
latter as a tube with a uniform circular section in the plane perpendicular to the tube centreline. 
The chain thickness, A, is defined as the radius of the circular section. Several frameworks 
have been introduced to capture the uniform thickness constraint in a way apt for numerical 
implementation. These approaches typically rely on a discretised representation of the thick 
chain 1 4-8 1491 1501 1281 15 II . In this study we shall employ the piece-wise linear modelling of 
the chain centerline introduced by Gonzalez and Maddocks |28|. 

We shall indicate with V the centerline of the chain consisting of a succession of points 
{r*o, ri, ...} equispaced at distance a. We shall further denote with hi the virtual bond joining 
the ith and i + 1th points, hi = r J+ i — ri. In order for the succession of points {ri} to 
be a viable centerline for a chain of thickness A, it is necessary that the radii of the 
circles going through any triples of points i, j and k, are not smaller than A. Accordingly, the 
Hamiltonian for the thick chain (tube) model can be written as 

n TC (r) = J2Mn jk ) (l) 

ijk 

where V3 is the three-body potential used to enforce the thickness A of the chain 1 28 49 52 
l29l l30l 1531 . As anticipated, the argument of V3 is the radius of the circle going through the 
triplet of distinct points i, j, k and has the form 

t/ / \ JO if r > A, 

W = { +oo otherwise. (2) 

Physically, the model of eqn.n m tr°duces conformational restrictions for the centerline 
that are both local and non-local in character, as depicted in Fig. [2 The local constraints 
are those where the triplet i, j and k identifies three consecutive points. The limitation on 
the radius of the associated circumcircle reflects the fact that, to avoid singularities, the local 
radius of curvature must not be smaller than A. This reflects on the following bound on the 
angle formed by two consecutive bonds: 

h^i±l > ! _ jL (3 ) 
a 2 ~ 2A 2 K ' 

On the other hand there is also a non-local effect due to the fact that any two portions of 
the centerline at a finite arclength separation cannot interpenetrate. It has been shown in ref. 
1 28 1 that this second effect can be accounted for by requiring that the minimum radius among 
circles going through any triplet of non-consecutive points, is also greater than (or equal to) 
A. The seamless way in which the local and non-local steric effects are accounted for make 
the model particularly appealing. Other discrete models relying on pairwise interactions for 
the excluded volume (such as the cylindrical model of refs. 1 44 45 1) may be adopted, though 
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ad hoc prescriptions for dealing with e.g. overlapping consecutive units need to be introduced 

Ell- 
in the present context we will consider the application of a stretching force, /, to the 
ends of a chain r = ro, rjv of thickness A (the contour length therefore being L c = Na). 
The Hamiltonian of eqn.^needs to be complemented with the stretching energy 

n = H v * ( r m )-f-(f N - ? ) ■ (4) 

ijk 

As customary we shall characterize the force dependence of the average normalised 
projection of the end-to-end distance, ?n — ro, along the direction of applied force: 

x = (5) 
Na\f\ 

where the brackets denote the canonical ensemble average. Owing to its self-avoiding nature, 
the stretching response of the chain cannot be characterised exactly by available analytical 
methods. We shall therefore resort to extensive Monte Carlo samplings, based on the 
Metropolis scheme, to evaluate the ensemble averages of eqn.[5] 

Besides the numerical study of the tube model subject to stretching it is interesting to 
illustrate a simplification of the model of eqn.^ which is amenable to an extensive analytical 
characterization. To do so we retain only the local thickness constraint and thus end up with 
a model that is essentially non-self-avoiding. The simplified nature of this problem, however, 
makes it very tractable also in the presence of a bending rigidity penalty, k^. We shall therefore 
consider the Hamiltonian 

ft = X! ^3(^+1^+2) - /• {r N - r ) - K b ^2 ~ t +1 > ( 6 ) 

i i 

again we stress that the three-body potential is restricted only to consecutive (local) triplets. 
In this form the stretching response of the model can be characterized exactly both at very 
low and very high forces using standard statistical-mechanical procedures 1 551 1561 . These 
two limiting regimes can be joined together to yield the following approximate expression for 
the stretching response of a locally-thick chain with bending rigidity (LTC+BR, for brevity): 



(3aF = 2K 




y(K, A/a) 



l+y(K,A/a) 2K ^ l + [l/2K) 

where K — Pki,, /3 — 1 /Kb T is the inverse Boltzmann factor and 

a 2 (_y 

2A 2 

coth(/-0 - ^ ' f < 0.5 



y(K,A/a)=i 1 &? {^* + *) f > °' 5 (8) 



(7) 



with z = -fgiK. The two cases in the above equation, reflect the fact that, for A < a/2 
no restriction applies to the angle formed by two consecutive virtual bonds since eqn. is 
always satisfied. 
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Expression[7]possesses some noteworthy limits. First, in the absence of thickness and in 
the continuum limit (a — > 0, K — > £ p /a, £ p being the persistence length) 1551 1561 . the model 
reduces to the well-known Marko and Siggia result for the WLC: 



1 



4(1 -x) 2 



(9) 



Secondly, in the absence of both thickness and bending rigidity one recovers the low- 
and high-force response of the freely-jointed chain with Kuhn length equal to a: 



3k B T , 



X 
X 



(10) 



It is of interest also the case of finite thickness A/a > 0.5 but no bending rigidity. In 
this case one obtains the following expression for the persistence length: 

a 



lnfl 



(11) 



4A 2 < 



Though this expression does not include the non-local self-avoiding effects it will be 
shown later to provide a good approximation of the persistence length obtained numerically 
for the full model of eqn.^ 

Local triplet 



Non-local 
triplet 




Figure 1. The finite thickness introduces steric constraints of local and non-local character 
that forbid configurations where the chain self-intersects. These constraints are conveniently 
treated within the three-body prescription of the thick-chain model. Within this approach the 
centerline of a viable configuration is such that the radii, R^jk of the circles going through any 
triplet of points on the curve i, j, k are not smaller than A. 

We conclude this section by mentioning that for the models of eqn.n] an d|6]the spacing 
of consecutive points is constant along the centerline, so that the contour length is unaffected 
by the application of forces of arbitrary strength. The inextensibility property is obviously a 
simplification of the behaviour found in naturally-occurring polymers which, at very high 
forces can undergo isomerization or structural transitions resulting in an "overstretching" 
beyond their nominal contour length. Several approximate treatments have been developed to 
correct the stretching response of inextensible models so to account for overstretching |[57ll5l 
by adopting additional parameters in the theory. In the present study we shall keep the number 
of model parameters to a minimum and hence postpone to future work the the investigation 
of the most suitable way to include overstretching in the TC model. 
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(a) (b) 

Figure 2. (a) Comparison between the persistence length of a thick chain obtained by Monte 
Carlo sampling (dotted-square) and that obtained from the local thickness approximation (solid 
curve) of eqn. <1 II . The inset illustrates the limitations of the local approximation for low 
values of A/a. (b) Tangent-tangent correlation (dotted-square) obtained from the analysis of 
10 4 (uncorrected) Monte-Carlo-generated chains of 1000 segments and with A/a = 2. The 
error bars which show the uncertainty of the plotted values are smaller than the square symbols. 
The solid curve is a single exponential fit to the numerical data yielding £ p = 16.71 ± 0.03. 

Numerical results 

The characterization of the stretching response of the thick chain of eqn. |4] was carried 
out using a Monte Carlo scheme: starting from an arbitrary initial chain configuration 
satisfying the thickness constraints, the exploration of the available structure space was done 
by distorting conformations by means of pivot and crankshaft moves. Newly-generated 
structures are accepted/rejected according to the standard metropolis criterion (the infinite 
strength of the three-body penalties of eqn.|2]was enforced by always rejecting configurations 
violating the circumradii constraints). 

The discretization length, a, was taken as the unit length in the problem and several 
values of A/a were considered, ranging from the minimum allowed value of 0.5 to the value 
of 4.0. This upper limit appears adequate in the present context since the largest nominal 
ratio for A/a among the biopolymers considered here is achieved for dsDNA for which one 
has A/a « 3.7 IH7I , For each explored value of A/a considered, we considered chains 
of length at least ten times bigger than the persistence length estimated through eqn. dl It . 
The relative elongation of the chain, was calculated for increasing values of the applied 
stretching force (typically about 100 distinct force values were considered). For each run, 
after equilibration, we measured the autocorrelation time and sampled a sufficient number 
of independent conformations to achieve a relative error of, at most, 10 -3 on the average 
chain elongation. For moderate or high forces this typically entailed the collection of 10 4 
independent structures while a tenfold increase of sampling was required at small forces due 
to the broad distribution of the end-to-end separation along the force direction. 

We first discuss the results for the persistence length obtained from the decay of the 
tangent-tangent correlations measured at zero force over an ensemble of sampled structures 
picked at times greater than the system autocorrelation time. The resulting data are shown in 
Fig. |2(a)| along with the curve corresponding to the approximate expression of eqn. ( fill . It 
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can be seen that the local approximation for the persistence length is very good in the range 
1.0 < A < 4, where the relative difference from the value found numerically is typically 
inferior to 10%. Significant relative discrepancies are, instead found as A approaches the 
limiting value of 0.5 (though it should be noted, the single exponential fit suffers from the 
very rapid decay of the tangent autocorrelation). In this case, only a narrow range of values 
for the angle formed by two consecutive bonds is forbidden. Consequently, the persistence 
length is very much affected by the (non-local) self-avoidance condition that is unaccounted 
for by the simple expression of eqn. dl \\ . As intuitively expected the value of £ p found 
numerically is larger than the one based on the local-thickness approximation. 
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(a) 



(b) 



Figure 3. Linear (a) and semi-log (b) plots of the force versus relative extension curves 
obtained from Monte Carlo sampling of chains of relative thickness A/a = 0.75, 1, 2 and 3. 



The analysis of the numerical results revealed several different stretching regimes in the 
elastic response of a thick chain. These are best discussed in connection with analogous 
regimes found in the two standard reference models, the freely-jointed chain (FJC) and 
the worm-like chain (WLC). We first point out that for both these models, as well as the 
generalisation of eqn. the relative elongation, x, depends linearly on the applied force, 
/. This result holds also for the thick chain model. However, due to the inclusion of the 
self-avoidance effect in the TC (absent in both the FJC and WLC), the Hookean relationship 
between / and x disappears in the limit of long polymer chains in favour of the Pincus regime, 



/ 



x 



i/(V"-i) 



3/5 being the critical exponent for self-avoiding polymers in three 



dimensions 1 25 16 58 1, see Fig.|4] 

For intermediate forces the Pincus behaviour is found to be followed by a second regime 
characterised by the same scaling relation found in the WLC at high forces, / ~ (1 — x)~ 2 . As 
shown in Fig. |4] at still higher forces the same scaling relation found in the FJC is observed, 
/ ~ (1 — a;) -1 . Physically, the first two regimes are determined by self-avoidance and chain 
stiffness or persistence length while the last regime is ascribable to the discrete nature of the 
chain l59ll6Ull55ll6Tl. 

In order to apply the thick-chain model to contexts where experimental data are available 
we have analysed the data of the numerical simulations with the purpose of extracting an 
analytical expression capturing the observed functional dependence of / on A and a. For 
any value of a and A the sought expression should reproduce the succession of the three 
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Figure 4. (a) Extension versus reduced force for a chain of relative thickness A/a = 1. Data 
points are presented in the (log(l — x), log(/)) plane to highlight the WLC- and FJC-like 
regimes found at moderate and high forces, (b) Illustration of the low-force crossover from 
the Hookean regime, x <x f, to the Pincus one, x <x / 2 / 3 for a chain of 1200 segments and 
A/a = 1. 



regimes discussed above. Among several trial formulae we found, a posteriori, that the best 
interpolation was provided by the following expression, 



/(*) 



ksT , / kxx 3/2 + k 2 x 2 + k 3 x 3 
tann 



a(l — x) 



1 



(12) 



where the parametric dependence on A and a is carried by the following expressions for the 
fc's, 



fcf 1 = -0.28394 + 0.76441 A/a 
0.15989 - 0.50503 A/a 
1.23330 A/a 



k 2 — 
k^ 1 = - 0.34984 



0.31858 A 2 /a 2 , 
0.20636 A 2 /a 2 , 
0.58697 A 2 /a 2 . 



(13) 
(14) 
(15) 



Within the explored ranges of A and a, the relative extension obtained from eq. H2i 
differs on average by 1% (and at most by 5%) from the true values at any value of the applied 
force, as shown in Fig. [5] 



Applications to experimental data and discussions 

We shall now discuss the application of the TC model to data sets obtained in DNA and 
RNA stretching experiments carried out for various [Na + ] concentrations. In particular, the 
data for dsDNA are taken from ref. Q (solutions of 10, 1.0 and 0.1 mM [Na + ]) while for 
single-stranded RNA (poly-U) we considered the recent results of ref. |62| (solutions of 
500, 300, 100, 50, 10 and 5 mM [Na+]). The fits of the experimental data is carried out 
through the standard procedure of minimizing the % 2 over the TC parameters: the contour 
length, L c , chain granularity, a and chain thickness, A. For the calculation of x 2 we nave 
estimated the effective uncertainty on the force measurements, taking the relative extension 
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(a) (b) 

Figure 5. Illustration of the performance of the parametric expression of eqn. 1121 in 
reproducing the stretching response obtained numerically for chains of relative thickness 
A/a = 0.75 and A/a = 4.0 (figures (a) and (b), respectively). In both cases the discrepancy 
between the computed and parametrised values is about 1% on average (and always smaller 
than 5%). 

as the independent variable, directly from the large data sets of poly-U (more than 400 points 
for each set). For dsDNA, owing to the more limited number of points (about 20-25) we 
propagated the declared experimental uncertainty on extension. 

The results of the fit procedures are provided in graphs[6]and0and tables^andlJ] 




5 10 15 20 25 30 35 

Extension (urn) 

Figure 6. Thick-chain fit of the experimental data (squares, circles and triangles) on dsDNA 
in solutions of various ionic strengths. The best-fit parameters are provided in TablelTI 

It is particularly appealing that, over the about 400 data points available for poly-U, the 
X 2 associated to the thick chain is very close to 1 for the set of measurements carried out for 
[Na + ] in the 50 to 500 mM range. In the case of the two smallest concentrations, [Na + ]= 
5 and 10 mM, a significant increase of the \ 2 is observed. The same is true for the fit of 
DNA measurements carried out in 0.1 mM [Na + ]. The worsening of the TC performance 
consequent to the increase of the screening length is reflected, among other effects, by the 
progressive importance of accounting for overstretching |63 6 62 1. Accordingly, the fitting 
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Figure 7. Application of the thick-chain model to Poly-U stretching data for (a) 300 mM 
[Na + ] and (b) 10 mM [Na + ]. Experimental data are shown as open circles, the fit with the TC 
is denoted with a solid line. The best-fit parameters for the TC are provided in Tablel2l 



dsDNA 


[Na+] 


a 


A 




L c 


x 1 


(mM) 


(nm) 


(nm) 


(nm) 


(nm) 




10 


9 ±2 


12 ±2 


55 ±2 


32.9 ± 0.2 


0.6 


1.0 


25 ±7 


26 ±4 


94 ±4 


32.3 ±0.2 


0.3 


0.1 


36 ± 16 


48 ± 13 


242 ± 14 


32.3 ±0.2 


1.9 



Table 1. Best-fit parameters obtained by applying the parametric force-extension expressions 
of the TC model to the experimental data on DNA in solutions of various ionic strengths. 



poly-U 


[Na+] 


a 


A 




L c 


x 2 


(mM) 


(nm) 


(nm) 


(nm) 


(nm) 




500 


1.10 ±0.01 


0.64 ± 0.01 


0.80 ±0.01 


2126 ±4 


1.49 


300 


1.09 ±0.02 


0.65 ±0.01 


0.91 ±0.01 


2123 ±5 


1.29 


100 


1.08 ±0.03 


0.68 ± 0.02 


1.11 ±0.01 


2117 ± 8 


1.05 


50 


1.11 ±0.03 


0.73 ± 0.02 


1.30 ±0.01 


2134 ± 8 


1.33 


10 


1.33 ±0.08 


0.94 ± 0.04 


1.94 ±0.03 


2138 ± 12 


4.21 


5 


1.41 ±0.10 


1.03 ±0.05 


2.26 ± 0.05 


2142 ± 12 


4.34 



Table 2. Best-fit parameters obtained by applying the parametric force-extension expressions 
of the TC model to the experimental data on poly-U in solutions of various ionic strengths. 

parameters obtained at the lowest salt concentrations can be expected to change upon the 
introduction of a stretching modulus. 

For both DNA and RNA the viability of the structural and elastic parameters of the TC 
fit can be compared against those obtained by different models and physical approaches. The 
most natural term of comparison for the elastic response is provided by the widely-employed 
WLC, which is the common reference model for determining the persistence length and 
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contour length of several types of biopolymers. It is appealing that both £ p and L c given 
by the TC for dsDNA and poly-U in Tables ^ and [2] are in close agreement with the same 
quantities given by the WLC, the relative difference being typically less than 5% for both 
quantities. For dsDNA, the values of £ p as a function of [Na + ] are compatible with the ones 
found by Baumann et al. (41 1 and those predicted by Lee and Thirumalai 1381 . The values 
of £ p for poly-U are, on the other hand, consistent with the WLC results found by Seol et 
al. in the original analysis of the data set provided to us |62|. Also, the order of magnitude 
of the granularity parameter, a, providing the best fit for various [Na + ] concentrations is 
consistent with the discretization length of the piece-wise cylindrical model optimised for 
DNA by Vologodskii and Frank-Kamenetskii 1641 . 

It is important to point out that the best fits of the TC model yield a \ 2 that is about 
one half the WLC one. This way reflect the fact that the WLC has one less parameter (the 
granularity) with respect to the TC model. However, the difference also stems from the very 
different functional dependence of f(x) in the two models. A noteworthy illustration of this 
fact is provided by the use of the Kratky-Porod (KP) model which represents a discretised 
chain with bending rigidity (e.g. a WLC endowed with a granularity parameter). The presence 
of the additional parameter, which allows the comparison with TC on a physically-equal 
footing, does not decrease appreciably the \ 2 which remains about twice the TC one. A 
considerable improvement over the WLC/KP fit is however possible upon the inclusion of 
the local thickness effects, that is through the LTC+BR model of eqns.[6]and0 The model is 
particularly interesting because it represents the simplest way of combining, in an approximate 
but analytical way, bending rigidity and thickness effects (within the piece-wise cylindrical 
model the inclusion of bending rigidity was recently discussed in ref. |45|). Despite the 
absence of non-local self-avoidance the LTC+BR model has a \ 2 performance that, on 
average, improves the TC one by 20%. In fact, for the three concentrations of [Na + ] reported 
in Table ^one obtains the following x 2 values: 0.6, 0.4 and 1.0. The associated thickness 
values (being approximately 14, 26, and 48 nm, for 10, 1.0 and 0.1 mM [Na + ], respectively) 
are consistent with the full TC determination, while the bending rigidity parameters are within 
50% of those of the "bare" WLC. These results indicate the benefit of supplementing the 
ordinary bending rigidity with the self-avoidance originating from the chain thickness. Owing 
to the fact that the combined model possesses a larger parameter space we will postpone to 
future work the detailed characterization of the model. 

We finally discuss the dependence of the effective thickness, A, on the concentration 
of counterions in solution (the apparent diameter of the polymer is simply given by 2A). 
For both DNA and poly-U, the effective radius shows a monotonic decrease for increasing 
concentrations and indicates reaching a saturated value (different for each type of the 
polymers), namely the inferred "bare" thickness. In particular, over the range of 0.1 to 10 
mM [Na + ] the apparent radius of dsDNA decreases from 48 nm to 12 nm. For ssRNA 
the decrease is instead, 0.96 nm to 0.63 nm over the range of 5 to 500 mM [Na + ]. As we 
anticipated in the introduction, the theoretical results derived by Stigter and Odijk 11341 1391 
have been previously used to predict the effective diameter of double-stranded DNA. The 
available theoretical predictions based on the approach of refs. 1341 l40l 1391 yield dsDNA 
radii of 8.1 nm for 10 mM [Na+], 27.1 nm for 1 mM [Na+] and 96.3 for 0.1 mM [Na+]. 
Given the very diverse nature of the approach and approximations of the present study and 
of refs. 1341 1391 . it is pleasing that the two determinations of the effective thickness are in 
reasonable agreement. 

We finally turn to the case of poly-U. Owing to its high flexibility, even short stretches of 
RNA cannot be modeled as stiff rods, and hence the theoretical analysis of ref. 1 34 1 cannot be 
used to generate a comparison term for the effective radius, A, of poly-U. It was pointed out by 
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Figure 8. Effective radius, A, of poly-U as a function of ionic strengths. The dashed curve is 
a visual guideline obtained from a spline interpolation of the points obtained from the TC fits. 

The approximate linear dependence of A on ■===. (with [Na+ ] expressed in molar units) 

V[iVo+] 

is shown in the inset. 



Stigter that for DNA the electrostatic contribution to the effective radius, though exceeding 
by several times the bare one, was of the same order of the Debye screening length, Ad 
(the "proportionality factor" ranging from about 4 to 8) |34|. Also in the case of RNA, 
the Debye length, which is proportional to 1/ ^/[Na+] is a useful concept for rationalising 
the dependence of the effective radius on the ionic strength. In fact, within the explored 
range of [Na + ], A appears to increase linearly with !./•»/ [Na+], as visible in the inset of 
Fig. |U(b). Assuming the validity of such linear relationship, the "bare" radius of poly-U 
is estimated by extrapolating A for vanishing Xd- One obtains Afc are w 0.59 ± O.Olnm, 
which compares well with the nominal value of about 0.5 nm observed in crystallographic 
structures of (non-hydrated) poly-U fragments (from PDB structure 1I5L). Incidentally we 
mention that also the value of a agrees with the nominal size of the modular poly-U units 
which is about 0.6 nm. Several efforts have been devoted in the past to establishing how the 
apparent persistence length of polyelectrolytes is affected by the interplay between the bare 
persistence length, the linear charge density and the ionic strength 1371 1361 l38l . Different 
manners of functional dependence of £ p on \r> are, in fact, observed for polyelectrolytes 
with different degree of flexibility, such as double-stranded and single-stranded DNA. It is, 
therefore, interesting to compare the dependence of the effective diameter of DNA and poly- 
U. The most notable difference is that the proportionality factor between the electrostatic 
contribution to the effective radius and the Debye length is of the order of 0.1, therefore 
appreciably smaller than for the case of dsDNA. 

Conclusions 

Extensive stochastic sampling techniques were used to characterize the behaviour of an 
inextensible thick polymer subject to a stretching force. The extension versus force response 
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found numerically was parametrised in terms of the polymer structural parameters (thickness 
and monomer length) and captured by an analytic expression. The resulting formula was 
used to fit experimental data obtained from stretching measurements of DNA and poly-U in 
solutions of different ionic strength. This represents a novel and physically-appealing route to 
extract the apparent structural parameters of polyelectrolytes starting merely from their elastic 
response. The inferred effective diameter for DNA was found to be in satisfactory agreement 
with the estimates obtained by Stigter through an unrelated approach. In particular, the 
electrostatic contribution to the effective DNA radius was found to be several times larger than 
the Debye screening length, Xrj. Also for the much-more flexible poly-U chain it is observed 
that the effective radius strongly depends on the ionic strength, having an approximately linear 
dependence on within the available ranges of [Na + ]. At variance with dsDNA, however, 
the electrostatic contribution to the effective radius was found to be an order of magnitude 
smaller than the Debye screening length. 
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